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SUMMARY 


A TAKOMAK Vacuum Vessel was analyzed using MSC/MASTRAN. Isoparametric quadri- 
lateral and triangular elements were used to represent the Vacuum Vessel shell 
structure. For toroidally symmetric loadings, MFCs were employed across model 
■boundaries and Rigid Format 2h was invoked, Un-symmetric loadings required the 
use of the cyclic symmetry analysis available with Rigid Format H 9 . NASTRAN 
served as an important analysis tool in the TOKAMAK design effort by providing 
a reliable means for assessing structural integrity. Interactive graphics were 
employed in the finite element model generation and in the post-processing of 
results. The authors -'eel that model generation and checkout with interactive 
graphics reduced the modelling effort and debugging man-hours significantly. 

INTRODUCTION 


This paper presents the structural modelling method, analysis procediores and 
results of a TOKAMAK Vacuum Vessel finite element analysis. The analysis 
described herein was required for the design-verification of the TOKAMAK 
FUSION TEST REACTOR (TFTR) Vacuum Vessel. The TOKAMAK is a toroidal device 
through which magnetic fields permeate to confine a plasma. The magneti'.' 
fields are produced by strong electric currents on the order of U0,000 amperes 
passing through both copper coils and the plasma. In the TOKAMAK, there are 
Gwo major magnetic fields: a toroidal field generated by current flowing in 

the coils enveloping the torus, and a poloidal field generated by current flow- 
ing through both the plasma and an equilibriinn field coil. By combining 
toroidal and poloidal magnetic fields, the TOKAMAK achieves a higher level of 
plasma stability than has been realized in any previous magnetic-confinement 
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systeii'. The increased plasma stability permits longer confinement times of 
higher temperatui'e plasmas. The TFTR will therefore come closer to meeting all 
necessary conditions for a net production of fusion energy than any previous 
magnetic fusion device. A basic function of the Vacuujii Vessel is to provide 
containlnent of the hot (100 million degrees Celsius) Deuterium-Tritium plasma 
while excluding the atmosphere. 


SYMBOLS 


Values ai'e given in both SI and US GustomEiry Units. Calculations were made in 
US Customary Units. 

fl 3 Bending stress, psi (MPa) 

imax Maximum combined stress, psi (MPa) 

fill Membrane stress, psi (MPa) 

fxy In-plEuie sheai’ stress, pal (MPa) 

Fx Circimiferential (toroidal) membrane load, Ib/in. (N/m) 

Fy Meridional membrane load, lb/ in. (N/m) 

Fxy In-plane membrane shear load, Ib/in. (W/m) 

Qx Transverse sheai’ load (causing Mx bending), Ib/in. (N/m) 

Mx Circiutiferential (toroidal) bending load, in Ib/in. (N-m/m) 

My Meridional bending load, in. Ib/in. (N-ra/m) 

qj. Transverse shear stress, psi (MPa) 

t Material Thickness,, in. (m) 

THE TOKAMAIC VACUUM VESSEL 

The Vacuum Vessel (Pig. l) is essentially a large doughnut-shaped vessel (or 
torus) composed of ten stainless steel segments joined at ten parting planes. 
Each segment contains five intersecting cylinders (three pie sections and two 
bellows cases) as shown in the p^lan view (Fig. 2). The diametral centerline 
of each cylindrical section intersects the vertical axis of syiranetry of the 
machine. Six of the segments include a Neutral Beam Injection Duct, as shoira in 
Figui'e 1. (The pui’pose of the Neutral Beam Injection Duct is to produce high- 
energy Deuterium atoms and to inject them into a magnetically-confined Tritiiuii 
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plasma causing Deuterium-Tritium fusion reactions to occur. ) Each segment is 
supported at both a radially inbofU'd and a radially outboard location. Inboard 
supports supply essentially vertical restraint, supporting approximately 1/3 of 
the deadweight of the Vacuum Vessel. Outboard supports allow free radially 
inward displacement while restraining radially outward, toroidal and vertical 
motion. 


FINITE ELEMENT IDEALIZATION 


Two models of Vacuum Vessel segments have been developed which span from part- 
ing plane to parting plane. A derivative of the non-Neutral Beam Injection 
Duct model was also created for analysis of anti-syinmetric loading. The first 
model, shown in Figure 3, contains detailed features such as diagnostic and 
access port extensions, non-structural port covers, the parting plane (P/P) 
weld thickness with its stiffening rings, the plasma limiter support structure, 
pie bellows-case (P/B) stiffening rings, and the inboard and outboard supports. 
QUADU and TBIA3 isoparametric plate bending elements were used for the shell 
idealization. (These elements were chosen because they provide superior 
results with fewer elements. The results of element test problems carried out 
by MSC (Ref. l) has shown that the acciiracy of the QUAD^) element is effectively 
independent of the aspect ratio. ) Beam elemeri,.s were used to simulate the 
stiffening rings at the P/B intersections, and those adjacent to the parting 
plane. These element types were also used to idealise the Vacuum Vessel sup- 
ports and flanges at the ends of the port extensions. Rod elements were used 
to model the protrusions at the port-shell intersections. Continuity of the 
finite element model across structural interfaces (viz. at the P/B intersec- 
tions, parting plane/ring intersections and at the vessel support connections) 
was achieved by means of rigid-bar and rigid-tri angular elements. The second 
model as shown in Figures 1* and 5 , is essentially the same as the first except 
for the inclusion of a Neutral Beam Injection Duct-. Each model consisted of 
app:.’oximately 1900 nodes with 6 DOF at each node. 
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MODEL GENERATION INTERACTIVE GRAPHICS 

The intex’active display featiu'e of Control Data' s UWISTRUC program facilitated 
the modelling of the numerous penetrations, eccentricities and extensions 
which comprise more than SQ% of the total model. In order to develop this 
geometrically complex finite element idealisation it was necessary to employ 
several pre-processors, (The sequence of steps employed in the model genera- 
tion which ore described below ai’e schematically depicted in Eigure 6.) First, 
a FORTRAN program was witten to geometi'ically outline the uoundai'ies of major 
individual sti'uctui'al components. The results of this program were input into 
UNISTRUC in teimis of "lines" and "points". The boundoi'y interiors were then 
meshed by employing UNISTRUC's generic element library (Ref. 2). 

In addition to nodes and elements , loads , material properties and physical 
properties were defined and graphically reviewed through UNISTRUC. This 
information was written to a "Neutral Input File" in an application-independent 
format. Inputting this file to the "NASTRAN Input-File Translator" resulted in 
a file containing pi’operly formatted Executive Control, Case Control and BuUc 
Data card images. To take advantage of more recent NASTRAN features which are 
not cxu'rently compatible with UNISTRUC (e.g. rigid elements) it was necessary 
to edit the cord images. This was accomplished through the use of the INTERCOM 
interactive editing featui'e of Control Data's Scope 3.*+ Operating System (Ref. 
3). Model bandwidth optiitiization aaid SEQGP card generation was accomplished by 
a stand-alone optimizer residing on a linked mainframe, and subsequently 
merged, through INTERCOM, with the NASTRAN BuHc Data. 

LOADINGS 

The TFTR is subjected to several different loading combinations. These com- 
binations include the Nomial Operating Condition, Bakeout, Discharge Cleaning, 
Plasma Disruption, Seismic and Coil Short-circuiting. As of this witing the 
Vacuum Vessel has been analyzed for the sevei-e conditions associated with a 
Plasma Disruption. A Plasma Disi'uption loading consists of electromagnetic 
centering forces, electromagnetic pressure and electromagnetically-induced 
toi'oidal and racking forces on the bellows case rings in addition to the 
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deadweight of the Vessel, and external pressurization. For the salce of 
brevity, this paper will only present results related to the gravity and 
pressurization subcase. 


AMALYSIS PROCEDURE 

Due to budgetary constraints it was decided early in the program not to model 
the entire torus structure. Instead, individual models consisted of a single 
segment type with boundary conditions selected so as to simulate a complete 
torus. Analysis with smaller models had shown that results of sufficient 
accuracy would be obtained with this idealization. For torcidally symmetric 
loadings, the most efficient analysis technique was to employ JiPCs across the 
boundaries of the selected segment type and invoke Rigid Format 24. Unsym- 
metric loadings required the use of the more expensive Cyclic Symmetry Analysis 
available in Rigid Format 49. (A Cyclic Symmetric run with a harmonic index of 
five was found to require twice the CPU time and 3.5 times the I/O as compared 
to the MPC'd boundary model.) Based on the economic and technical resources 
available, it was decided to perform the analyses on a second mainframe con- 
taining the MSC/NASTRAN Version 40 operating under the Network Operating 
System/Batch Environment. 


ANALYSIS RESULTS 


The resulting significant internal membrane and bending force distributions are 
depicted in Figures T through 10. Continuity of the curves was obtained by 
smoothing the computer-output values at element centroids . In regions of load 
and/or structirral discontinuity, the results were extrapolated to the edges of 
the finite elements vmder consideration. These results were checked, where 
feasible, by utilizing equilibrium and compatibility relationships, thus 
maximizing the information yielded by the analysis. Post processing of the 
results would have been mere expeditious had an option been available in 
NASTRAN to print out element internal load intensities (Mx, My, Nx, Ny, Nxy, 

Qx, Qy) at the element corners in addition to the element centroid. This 
feature would more clearly define internal load gradients. Had this option been 
available , much of the time expended in manually smoothing and extrapolating the 

369 



results would have been saved. Figures 7 and 8 present the critical internal 
loads below the torus' horizontal :;.xis of syinmetry, specifically at the follow- 
ing critical locations; the raaocimum pie /bellows interface eccentricities and 
at the connection of the outboard supports to the torus , Figures 9 and 10 each 
present the internal loads along a meridian. Figure 9 depicts the axial and 
bending loads along the centerline of pie 2, while Figure 10 depicts the loads 
in the parting plane weld. 

The relationship between the internal loads and the resulting stresses is as 
follows : 
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combined membrane and bending 
stresses 


In-plane shear stress 





Transverse shear stress 


where t = .5 inches (.0127 meters) applicaWe to Figures 7 thru 9 

t = .312 inches (.00792 meters) applicable to Figure 10, parting 
plEuae weld 


An inspection of the results indicates that the critical stresses occur ir the 
following Vacuum Vessel locations: 


-In Pies 1 and 2, at the maximum pie/bellows interface eccentricities. 
Figure 7. 


- Ik * = -9Q0 „ 6 X 450 . ^ -12,600 psi 

^ max " t %2 .5 ^ ^j2 (-86.874 MPa) 

-In Pie 2, at the diagnostic port-shell intersection. Section A-A, 
Figure 8, 
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6m 

X 



-4500 

.5 


6 X 100 _ -11,400 psi 

5“ “ (-78.6 MPa) 

(. 5 )^ 


370 



- In Pie 2, at the diagnostic port-shell intersection, Figure 9, 


fy 


max 



-2500 

.5 


6 X 130 

(. 5 )^ 


-8120 psi 
(-55.9 MPa) 


In the parting plane weld 
creased hy 20S? to account 
inch (.0079^ meters) weld 


fx 


max 



1.2 


region, Figure 10, the Mx bending results are in- 
fer moment peaking at the intersection of the 5/18 
with the p/p stiffening ring: 

^^x -750 ^ ^ 6 X 100 _ -9800 psi 

" .312 " ■ TSI? ” (-67.57 MPa) 


DEFLECTIONS 


Grid point displacements, element forces and stress were output to the NASTRM 
UTl file by appropriate DMAP instructions generated by UNIST3UC. The contents 
of the UTl file were transformed by means of the UNISTRUC file translator for 
compatibility wit’. e UNISTRUC system, thus making interactive display of the 
results possible. Displays of the deflected shapes of major structural 
components appear in Figures 11 thru l6. A maximum deflection of .037 inches 
(0.0009^i meters) oc cured at the intersection of the outside edge of the 
diagnostic port and the shell on Pie 2 as shown in Figure l4. 


CONCLUSIONS 


NASTRAN proved to be a valuable analysis tool for the design-verification of a 
TFTR Vacuum Vessel. Furthermore, it was found that analysis of a 1/10 segment 
of a structure with Rotational Cyclic Symmetry (K = 5) and non-symmetrical 
loading was approximately 2 1/2 times more costly than for a syinmetrically 
loaded model of the same fundamental region with MFCs across the boundaries. 

The use of interactive graphics in hoth the pre and post-processing modes was 
found to he an indispensible tool when dealing with complex three-dimensional 
models . 


371 



REFERENCES 


1. MacNeal-Schwendler Corporation, MSC/MASTRM Application Manual For 
CPC 6000 Series , Los Angeles, California. 

2. Control Data Corporation, Unified Structural Design System (UNISTRUC) , 
Version 1.2 , Publication No. T6079600, Publications and Graphics 
Division, St. Pa\il, Minnesota. 

3. Control Data Corporation, CYBERNET Service INTERCOM U Reference Manual , 
Publication No. 8 U 0 OOO 2 U, Publications and Graphics Division, 

Data Services Publications, Minneapolis, Minnesota, 


372 



ORIGINAL 
OF POOR QUALITY 



tb) Vacuum Vessel 


2036-001 


Figure 1 Tokamak Fusion Test Reactor - Vacuum Vessel 
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Figure 3 TFTR Vacuum Vessel - Finite Element Model, Isometric View 
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'Figure 4 Neutral Beam Duct Segment - Ftnite Element Model, Plan View 



Figure 5 Neutral Beam Duct Segment * Finite Element Model, Isometric View 
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Figure 6 Pre and Post Processing with Interactive Graphics 
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C,PP STIFF RING 




10'" LB/IN (175 X 10*’ M/M) 



Figure 9 Loads Along the Centedine of Pie 2 
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LOAD CONDITION; 

PRESSURE & GRAVITY 

'"shell = -25 PSI (-.1 724 MPa) 

'"bellows case = + 5 PSI (.0345 MPa) 
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Figure 11 Deformed Sh?ipe of TFTR Vacuum Vessel Finite Element Model, Isometric View 
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Figure 13 Deformed and Undeformed Profiles of Bellows Case 2 
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